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Abstract 


A  model  developed  previously  to  investigate  the  transport  of 
electromagnetic  fields  and  energy  in  solid-armature  raUgims  is 
simplified  to  exclusively  study  heating  within  the  interior  of  the 
armature.  A  particular  t5^e  of  armature  is  then  specifically  studied, 
and  the  effect  of  variation  of  a  number  of  parameters  that  affect  its 
geometry  is  investigated.  Results  of  the  calculations  are  studied 
carefully  and  explained  by  basic  physical  principles.  Some  indication 
of  how  the  results  and  conclusions  can  be  used  to  influence  design 
criteria  for  armatures  is  also  indicated. 
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CURRENT  AND  HEAT  TRANSPORT  IN  A 
DOUBLE  TAPER  SABOT-ARMATURE 


1.  Introduction 


In  previous  work  (Powell,  Walbert,  &  Zielinski  1993),  we  developed  a  two- 
dimensional  (2-D)  model  for  investigating  the  coupled  effects  of  electromagnetic 
field  transport  and  heat  transport  in  solid-armature  railguns.  The  model  has  been 
extended  on  several  occasions  in  recent  years,  and  studies  have  been  undertaken 
for  various  types  of  armatures  (Powell  &  Zielinski  1995a,  1995b,  1997,  1999). 
These  studies  have  enhanced  our  understanding  of  how  force  and  thermal 
loading  depend  on  geometry,  how  velocity  affects  the  manner  in  which  current  is 
conducted,  and  how  basic  physical  phenomena  can  be  used  to  influence  design 
criteria.  More  recently,  results  from  the  model  calculations  have  been  employed 
in  structural  dynamics  codes  to  predict  the  response  of  both  the  armature  and  the 
rails  to  the  effects  of  electromagnetic  acceleration  (Newill,  Zielinski,  &  Powell 
2001). 


A  schematic  diagram  of  a  cross  section  of  a  railgvm,  which  has  a  simple 
rectangular  armature,  is  shown  in  Figure  1.  The  two  dimensions  shown  are  the 
actual  dimensions  treated  in  the  calculation.  Current,  conducted  in  the  direction 
indicated  by  the  arrows,  produces  a  magnetic  induction  field  B  (which  points 
out  of  the  page)  in  the  inner  bore  of  the  gun.  This  induction  field  interacts  with 
the  current  in  the  armature  via  the  standard  Lorentz  force  to  accelerate  the 
projectile  down  the  barrel  of  the  gun.  Current  initially  flows  along  the  inside 
surfaces  of  both  the  rails  and  the  armature,  but  it  eventually  diffuses  into  the 
interior  of  those  conductors.  The  intent  of  the  model  calculations  is  to  determine 
the  temporal  and  spatial  distribution  of  the  electromagnetic  fields,  current 
density,  and  energy  within  the  armature  and  rails. 


Models  similar  to  the  one  discussed  here  have  been  developed  in  the  past  and 
summarized  extensively  in  earlier  work.  The  most  closely  related  work  was 
probably  the  2-D  analysis  undertaken  by  Long  with  both  analytical  (1986)  and 
finite  element  (1987)  models.  Those  models  did  much  to  expand  our 
understanding  of  how  various  parameters  affected  the  operation  of  the  railgvm. 
Later,  finite  element  models  were  developed  in  three  dimensions  by  Rodger, 
Leonard,  and  Eastman  (1991)  and  by  Hsieh  (1995).  The  model  by  Hsieh,  known 
as  EMAP3D,  has  been  vmdergoing  particularly  active  development  during  the 
last  decade  approximately  and  offers  promise  for  explaining  many  outstanding 
problems  in  railgvm  dynamics. 
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Rail 


Armature  Projectile 


Figure  1.  Schematic  Diagram  of  Railgim. 


In  the  present  work,  we  develop  a  simplification  of  our  previous  models,  which 
will  be  used  exclusively  to  examine  current  and  heat  transport  in  the  armature. 
The  rails  are  not  considered  in  the  calculations.  Consequently,  we  do  not  accoimt 
for  important  phenomena,  such  as  the  velocity  skin  effect,  that  occur  at  the  rail- 
armature  interface.  While  accounting  for  these  effects  is  clearly  necessary  in 
modeling  the  complete  launch  problem,  neglecting  them  is  probably  satisfactory, 
provided  we  confine  our  attention  to  a  region  that  is  well  within  the  interior  part 
of  the  armature.  The  rationale  for  adopting  this  simplification  is  that  it  provides  a 
mechanism  by  which  calculations  can  be  done  in  a  reasonably  simple  way  and  in 
a  reasonably  short  time. 

We  will  furthermore  confine  our  attention  to  a  type  of  armature  known  to  be  of 
great  value  in  laimching  high  der\sity  tactical  payloads.  An  illustration  of  the 
cross  section  of  such  an  armature,  derived  from  engineering  criteria,  is  shown  in 
Figure  2a  (Zielinski  2001);  a  more  detailed  but  schematic  illustration  of  what  is 
actually  modeled  is  shown  in  Figure  2b.  The  armature  has  reflection  S5mmetry 
about  the  centerline  at  y  =  0  so  that  it  is  necessary  to  consider  only  the  upper 
half.  The  left-hand  face  is  formed  by  two  lines,  L,  and  L^,  which  make  angles  a 
and  P  with  the  horizontal,  and  by  an  arc  of  radius  r„.  The  center  of  the  arc  is  at  the 
point  {x^}/g),  and  lines  Lj  and  are  tangent  to  the  arc  at  the  points  and 

(x^y^).  The  upper  left-hand  comer  of  the  armature  is  taken,  for  convenience,  to  be 
at  x  =  Q ,  whereas  the  lower  left-hand  comer  is  at  x  =  Xy  Current  is  conducted 
across  the  centerline  at  y  =  0,  throughout  the  interior  of  the  armature,  and  across 
the  top  at  y  =  h.  We  are  interested  in  determining  how  a  variation  in  certain 
dimensional  parameters  affects  the  spatial  and  temporal  distribution  of 
electromagnetic  fields,  electromagnetic  forces,  and  energy  within  the  interior  of 
the  armature. 


a.  Illustration  of  cross  section. 


(0,h) 


b.  Schematic  depicting  details  of  the  model. 
Figure  2.  Armature  Employed  in  the  Analysis. 


The  organization  of  the  report  is  as  follows.  In  Section  1,  we  describe  the  model 
in  some  detail  and  derive  the  governing  equations  that  need  to  be  solved.  In 
Section  3,  we  present  and  compare  results  of  some  specific  calculations  in  which 
the  input  parameters  a,  p,  and  are  varied  in  a  range  of  values.  In  Section  4,  we 
present  conclusions  and  indicate  specific  problems  that  require  further  study. 


2.  Model  and  Governing  Equations 


As  indicated  previously,  we  employ  a  simplification  of  the  model  that  we  have 
used  before  to  study  current  and  heat  transport  in  the  armature  and  rails  for 
solid-armature  railgims.  The  model  is  two  dimensional  in  the  x-y  plane,  as 
shown  in  Figure  1.  We  are  interested  in  determining  the  electromagnetic  fields  as 
well  as  the  temperature  within  the  armature  as  a  function  of  both  space  and  time. 
The  armature  carries  a  total  current  per  unit  height  (i.e.,  distance  along  the  rail  in 
the  z  direction)  given  by  j( t).  That  current  enters  at  the  bottom  at  y  =  0  and 
leaves  at  the  top  at  y  =  h. 

2.1  Governing  Equations  and  Boundary  Conditions 

As  in  previous  work,  it  is  convenient  to  work  in  a  frame  of  reference  in  which  the 
armature  is  at  rest.  The  basic  equations  that  determine  the  electromagnetic  fields 
and  temperature  consist  of  Maxwell's  equations  as  well  as  an  energy 
conservation  equation.  Let  B ,E ,  and  J  be  the  magnetic  induction,  the  electric 
field  intensity,  and  the  current  density,  respectively.  Let  p  and  a  be  the  magnetic 
permeability  and  the  electrical  conductivity,  respectively.  Then,  from  Maxwell's 
equations,  we  have 


VxE  =  -^,  (1) 

Bt 

(2) 

and 

J=gE.  (3) 

Equations  (1)  and  (2)  are  the  standard  Maxwell  curl  equations,  and  Equation  (3) 
is  Ohm's  law.  In  writing  Equation  (2),  we  have  neglected  the  displacement 
current.  In  general,  this  approximation  is  valid  for  the  relatively  large  time  scales 
in  which  these  types  of  conduction  problems  occur. 

In  accordance  with  the  symmetry  conditions  as  well  as  the  2-D  nature  of  the 
model,  we  assume  that  =  0,  E^=0 ,  B  =  Bd^  and  that  none  of  the  variables 
depend  on  the  z  coordinate.  With  these  assumptions,  it  is  possible  to  uncouple 
Equations  (1)  through  (3)  to  produce  a  single,  second  order,  differential  equation 
that  governs  the  diffusion  of  the  induction  field  B.  The  result  is 

a-g  1  da  dr  dB  1  dcr  dr  dB 

^  dt  dx^  dy^  G  BT  Bx  Bx  G  BT  By  By 
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In  obtaining  the  result,  we  have  assumed  that  a  depends  on  position  only 
through  the  temperature  T. 

To  accoimt  for  energy  transport  and  ohmic  heating,  we  assume  for  the  energy 
conservation  equation  that 


^dT  d^T  d^T  dK 

pC—=K—T+^^rT-  +  ^ 
dt  dx^  dy^  dT 


rart  d/cfdr)'  i 


V  y 


ar 


dy 


rasY 

-  +  - 

dx  J  ay 


-Bl-  (5) 


In  this  equation,  p,  C,  and  ic  represent  the  density,  specific  heat,  and  thermal 
conductivity,  respectively,  and  these  variables  are  also  assumed  to  be  functions 
only  of  the  temperature.  The  terms  involving  the  induction  field  in  the  equation 
can  be  recognized  as  just  7^/(7  and  therefore  correspond  to  the  energy 
dissipated  per  unit  volume  and  time  as  a  result  of  ohmic  heating.  The  term 
represents  the  energy  absorbed  during  melting  and  depends  on  the  heat  of 
fusion  of  the  armature  material.  A  specific  form  for  this  term  is  discussed 
subsequently.  In  writing  the  equation  in  this  form,  we  have  assumed  that  ohmic 
heating  is  the  only  source  of  energy  generation  and  that  thermal  conduction  is 
the  only  energy  transport  mechanism.  As  implied,  melting  is  accoimted  for,  but 
no  other  phase  change  is  assumed  to  occur.  It  is  evident  that  Equations  (4)  and 
(5)  are  coupled  since  the  electrical  conductivity  depends  on  temperature,  and  the 
source  term  in  the  energy  equation  depends  on  the  induction  field. 


Equations  (4)  and  (5)  are  identical  to  those  employed  in  our  previous  work, 
except  that  there  are  no  terms  that  account  for  the  relative  motion  of  the 
armature  and  rails.  The  absence  of  any  velocity  effects  is  because  the  rails  are  not 
included  in  the  calculation  and  because  the  equations  are  solved  in  the  moving 
reference  frame.  The  equations  must  be  solved  subject  to  bormdary  conditions  on 
all  sides  of  the  armature.  We  assume  that  there  is  no  transfer  of  heat  from  the 
armature  to  its  surroundings  so  that  at  every  surface  we  must  have 


n»Vr  =  0,  (6) 

where  n  represents  the  unit  normal  to  the  surface  in  question.  On  the  left-hand 
face  of  the  armature,  we  assume  that 


B  =  M/,  (7) 

and  on  the  right-hand  face  that 

5  =  0.  (8) 

Along  both  the  lower  and  upper  surfaces,  we  assume  that 
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With  the  exception  of  the  conditions  along  the  upper  siurface,  these  boundary 
conditions  are  similar  to  those  that  we  have  justified  previously  and  employed  in 
all  our  calculations.  The  conditions  on  the  upper  face,  where  in  reality  the 
armature  is  in  contact  with  the  rail,  are  not  realistic.  However,  since  we  are 
primarily  interested  in  analyzing  heating  at  points  well  within  the  interior  of  the 
armature,  we  assume  that  file  precise  conditions  at  that  interface  do  not  greatly 
affect  the  results.  For  simplicity,  we  have  chosen  the  conditions  to  be  the  same  as 
those  along  the  centerline  at  y  =  0. 

2.2  Coordinate  Transformation 

As  in  previous  work,  we  find  it  convenient  to  transform  the  equations  of  the 
previous  section  to  a  coordinate  system  in  which  the  space  occupied  by  the 
armature  is  rectangular.  Various  transformations  with  various  degrees  of 
complexity  are  possible,  but  for  this  problem,  we  choose  a  simple  algebraic 
transformation  to  new  coordinates  ^  and  t|  given  by 

^  =— 

L-s 

71  =  y  (10) 

T  =  t 


where  s  =  s{y)  represents  the  horizontal  distance  from  A:  =  0to  the  left-hand 
edge  of  the  armature  at  some  point  y.  It  is  evident  that  this  transformation  will 
map  the  left-hand  edge  of  the  armature  into  <|  =  0,  while  the  right-hand  edge  is 
maintained  at  ^  =  L.  For  the  particular  type  of  armature  being  considered,  s  is 
given  explicitly  by  the  relations 

(JC,  -Xj)  f 

s  =  — - fory<y2 

^2 

^  =  +  for  >^2  ^  ^  >1 

h-y  c 

s  =  - — fory>yi. 

h-y, 


The  coordinates  Xj,  y^,  and  so  forth  denote  coordinates  of  the  points  indicated 
in  Figure  2b.  In  general,  we  only  specify  the  parameters  a,  h,  y^  and  L  as 
input.  The  remaining  parameters  follow,  however,  from  a  geometrical  analysis  of 
Figure  2b;  formulas  for  these  quantities  are  given  in  Table  1, 


Table  1.  Coordinates  of  Points  Indicated  in  Figure  1 


Parameter  Formula  in  Terms  of  Input  Parameters 


Xq  (h- yQ-rQCOsa)/tana-rQsma 

x^  (/z-Jq -r(,cosa)/tana 

Ti  Jo  +  rocosa 

X,  {h  -  yQ-r^  cos  a) /tan  a -r^  sina:  +  sin  P 

T2  To-'bCOSjS 

(h  -  -  To  cos  a) /tan  a -Tq  sinQ:  +  r,,  sin  P  -(y^-  cos  /3)/tan 


Equations  (4)  and  (5)  are  transformed  to  the  new  coordinates  by  the  application 
of  the  chain  rule  to  express  derivatives  with  respect  to  x,  y,  and  t  in  terms  of 
derivatives  with  respect  to  r\,  and  x.  After  considerable  algebra,  we  find  that 
the  equations  become 


H<7  —  =  A 


2  , 


dr 


■+r 


ds 

y9yj 


d-B  d^B  ^  ds  d^B  d^sdB 
d^^^dri^  ^dyd^dri  ^  dy^  d^ 


Ixif  da  dT  dB  do  f  ds 


L-s 


yiy, 


d^-adTd^d^  cjarlayla^a^  odTdridr] 


dT  dB  1  da  dT  dB  (12) 
■  + 


y/  do  ds  dT  dB  y/  dcr  ds  dT  dB 


and 


^ dT  aio  1  a  1  -  os  a  1  - 

pC—  =  KA^^+K^-2YK—-—+Ky/^ 

dr  d^"  dr}"  dy  d^dr] 


d^T  d^T  .  ds  d^T  .fds'^' 


dy 


d^T  d-s  dT 
2-W - 


Kyr 


L-s 


j 


^  ar 


2  dK 

+y/  — 
^  dT 


dT 

a<f 


— 


a/  a^ 

dr]  j 


2  A- 


(13) 


dT  dy  dr]  d^  fX^G 
ly/  ds  dB  dB 


fas'! 

2  2 

+  ^ 

^asY 

fas] 

^  1 

A 

f— 1 

1 

|ro 

IT  a 

[ayj 

i^lj 

*  2 
ji  a 

,ai7, 

[x^G  dy  dr]  d^ 

In  these  equations,  A  and  v|/  are  given  by  the  relations 
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and 


¥  =  k-^-  (15) 

L-s 

Similar  transformations  must  be  applied  to  the  boundary  conditions,  but  those 
derivations  are  straightforward  to  perform. 

2.3  Material  Properties  and  Latent  Heat 

In  all  the  calculations  undertaken,  the  armature  was  assumed  to  be  composed  of 
aluminum  (7075).  Material  properties  were  obtained  from  curve  fits  discussed 
previously  (Powell  &  Zielinski  1997).  In  those  fits,  it  was  assumed  that  the 
specific  heat  C,  the  resistivity  77^,  the  thermal  conductivity  k,  and  the  density  p 
could  be  represented  as  linear  functions  of  the  temperature  T.  Consequently,  for 
any  function  F,  we  have 


F  =<3o+a,r.  (16) 

Values  of  and  a,  for  solid  and  liquid  Al(7075)  are  shown  in  Table  2.  All  values 
are  in  the  appropriate  international  system  (SI)  units. 


Melting  was  assumed  to  occur  between  values  of  temperature  given  by 
7].  =  750  K  and  =  908  K.  At  temperatures  below  T,,  solid  properties  were 
used,  while  at  temperatures  above  T^,  liquid  properties  were  used.  At 
intermediate  values  of  T,  i.e.,  for  the  solid-liquid  mixture,  the  property  in 
question  was  scaled  linearly  with  T  between  the  value  for  the  solid  at  T,  and  for 
the  liquid  at 

Table  2.  Values  of  Coefficients  for  Material  Properties  of  Al(7075) 


Property 

Solid 

^0 

Liquid 

c 

645.0 

0.705 

1302.0 

0 

Vr 

1.42x10'® 

1.36x10'’“ 

1.51x10'' 

1.33x10'’“ 

K 

247.0 

-6.92x10'' 

93.0 

0 

P 

2851.0 

-0.5 

2294.0 

0 

During  the  time  that  melting  is  occurring  at  some  location  in  the  armature,  an 
amount  of  energy  per  unit  volume  and  time  is  required  to  convert  solid  to 
liquid.  Generally,  this  energy  is  represented  by  the  relation 
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E,=pLF(T)^ 

at 


(17) 
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where  L  denotes  the  latent  heat  of  fusion  and  F  is  a  function  that  corresponds  to 
the  amount  of  liquid  produced  per  unit  rise  in  temperature  of  the  solid-liquid 
mixture.  For  most  metallic  elements,  melting  takes  place  at  a  well-defined 
temperature  so  that  F  can  be  represented  as  a  delta  function,  5{T  -T^),  in  which 
is  the  melting  temperature.  For  most  armature  materials,  however,  melting 
occurs  over  a  range  in  temperature,  i.e.,  between  T.  and  T^,  and  the  function  F  is 
not  known  with  any  accuracy.  We  make  the  same  assumption  as  in  our  previous 
work  that  this  function  can  be  represented  as 


F{T)  = 


1  £ _ 

niT-TJ^+e^^ 


(18) 


where  =  (T-  +Ty)/2  and  e  is  chosen  so  that  the  function  suitably  peaks  about 
T^,,.  This  particular  choice  is  not  unique,  but  it  seems  reasonable  until  more 
definitive  information  is  known;  furthermore,  it  approaches  the  delta  function  in 
the  limit  that  e  approaches  zero.  In  our  actual  calculations,  we  chose 
£  =  {T^  —  7))/6;  we  have  demonstrated  previously  that  this  choice  implies  that 
about  90%  of  the  latent  heat  is  absorbed  in  the  relevant  temperature  range.  The 
latent  heat  is  taken  as  £  =  4x10^  J  /kg. 

2.4  Numerical  Solution 

The  numerical  solution  of  Equations  (12)  and  (13)  is  performed  in  the  following 
manner.  A  non-uniform  rectangular  grid  in  the  computational  space  of  ^  and  Tj  is 
generated.  The  grid  points  are  clustered  so  as  to  achieve  a  relatively  fine  grid  in 
the  physical  space  of  x  an  y  in  regions  where  significant  spatial  variations  occur. 
Such  a  region  might  correspond,  for  example,  to  the  left-hand  edge  of  the 
armature  where  diffusion  is  initiated  and  where  substantial  curvature  in  the 
surface  occurs.  The  equations  are  then  solved  in  the  computational  space,  and 
the  results  are  again  transformed  into  the  physical  space  for  analysis  and 
presentation. 


The  numerical  technique  that  is  applied  to  solve  Equations  (12)  and  (13)  has  been 
described  in  earlier  work  (Powell,  Walbert,  &  Zielinski  1993)  and  is  not  described 
in  detail  here.  In  brief,  all  derivatives  are  represented  by  standard  finite 
differences,  and  the  resulting  nonlinear  equations  are  solved  by  iteration  (Ames 
1997).  Some  specific  grids  employed  are  discussed  subsequently.  Typically,  a 
time  step  of  50  to  100  ns  proves  satisfactory. 
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3.  Calculations 


We  now  describe  results  of  seven  calculations  in  which  we  have  varied  the  input 
parameters  a,  P,  and  and  have  endeavored  to  determine  the  effect  of  that 
variation  on  heating  in  the  vicinity  of  the  root  radius  [i.e.,  to  In  all 

calculations,  the  other  input  parameters  were  kept  constant  at  the  following 
values:  L  =  10cm,  A  =  38  mm,  and  yo  =  18mm.  Values  of  a,  P,  and  are 
summarized  for  the  various  cases  in  Table  3.  In  Cases  2  and  3,  was  varied 
relative  to  the  value  in  the  baseline  case.  Case  1,  while  a  and  p  were  held 
constant;  in  Cases  4  and  5,  the  angle  a  was  varied  while  and  P  were  held 
constant;  and  in  Cases  6  and  7,  the  angle  P  was  varied  while  and  a  were  held 
constant.  Also  indicated  in  the  table  are  values  for  the  other  parameters  in 
Table  1  (i.e.,  x,,  x^,  y^,  y^,  Xydind  x^  that  change  as  a  result  of  the  variation  in  a,  P, 
and  Tg 


Table  3.  Input  Parameters  for  Various  Calculations 


Input 


Parameter 

Case  1 

Case  2 

Case  3 

Case  4 

Case  5 

Case  6 

Case  7 

Tg  (mm) 

3.5 

2.0 

5.0 

3.5 

3.5 

3.5 

3.5 

a  (deg.) 

18 

18 

18 

12 

24 

18 

18 

P  (deg.) 

12 

12 

12 

12 

12 

8 

16 

Xj  (mm) 

51.31 

55.70 

46.92 

77.98 

37.74 

51.31 

51.31 

x^  (mm) 

50.96 

55.50 

46.41 

77.98 

37.04 

50.71 

51.19 

y,  (mm) 

21.33 

19.90 

22.76 

21.42 

21.20 

21.33 

21.33 

yjmm) 

14.58 

16.04 

13.11 

14.58 

14.58 

14.53 

14.64 

Xg(mm) 

50.23 

55.08 

45.37 

77.26 

36.32 

50.23 

50.23 

Xj  (mm) 

-17.62 

-19.98 

-15.26 

9.41 

-31.53 

-52.70 

0.15 

The  current  per  unit  height  along  the  z  direction  (i.e.,  the  height  of  the  rail)  was 
also  the  same  for  each  calculation  and  was  taken  to  be  of  the  form 

j=  joSinint/to).  (19) 

In  this  equation,  tg  is  the  acceleration  time  and  the  time  when  the  current  goes  to 
zero.  The  amplitude  jg  was  chosen  so  that  this  2-D  model  would  have  produced  a 
final  velocity  of  the  armature  of  for  a  given  h  and  m',  where  m'is  the  armature 
mass  per  unit  length  along  the  z  direction.  For  the  assumed  values  of  A  =  38  mm, 
m'  =  20  kg  I  m,  and  =  2.3  km  /  s,  we  found  that  jg  =  2.7  x  1 0’  A  /  m. 
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Reasonable  corroboration  between  modeling  and  experimental  results  has 
previously  been  obtained  via  this  scaling  procedure  (Powell  &  Zielinski  1995). 
For  a  real  3-D  railgun,  in  which  the  in-bore  fields  are  smaller  on  average  than  for 
the  2-D  case,  a  higher  value  of  (by  about  a  factor  of  2)  would  be  required  to 
achieve  the  same  acceleration. 

Shown  in  Figure  3  is  a  t)q)ical  physical  space  grid  employed  in  the  calculation  for 
Case  1.  As  was  pointed  out  previously,  this  grid  was  generated  from  a 
rectangular  computational  space  grid  that  had  a  non-uniform  grid  spacing.  In 
particular,  the  grid  was  chosen  so  that  was  small  near  the  left-hand  edge  of 
the  armature  where  diffusion  initiates,  and  Ar\  was  small  in  the  vicinity  of  the 
root  radius.  The  grid  shown  in  Figure  3  contains  80  points  in  the  x  direction  and 
70  in  the  y  direction.  We  have,  however,  used  both  larger  and  smaller  grids  to 
perform  the  calculations.  The  largest  grids  employed  contained  about  2,500  grid 
points  and  gave  satisfactory  results  except  at  very  small,  localized  regions  where 
there  was  a  transition  from  a  straight  line  segment  to  a  curved  segment.  These 
discrepancies  could  probably  have  been  eliminated  with  fewer  points  clustered 
more  judiciously. 


Figure  3.  Physical  Space  Grid  for  Case  1. 


In  Figure  4  are  plotted  lines  of  constant  magnetic  induction  at  the  time  of 
maximum  current,  namely,  1.3  ms.  It  is  possible  to  prove  from  Maxwell's 
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equations  that  current  cannot  cross  these  lines,  so  they  can  also  be  taken  to 
represent  current  streamlines.  Eleven  such  streamlines  are  shown  in  the  figure 
(counting  the  left-  and  right-hand  boundaries),  so  10%  of  the  total  current  is 
contained  between  successive  lines.  It  is  evident  that  the  current  has  not  yet 
diffused  through  the  armature  at  this  early  time  and  is  still  confined  largely  to 
the  left-hand  side.  The  current  density  is  high  in  the  vicinity  of  the  root  radius,  as 
can  be  seen  from  the  proximity  of  the  streamlines.  This  high  current  density 
results  from  current  being  conducted  aroxmd  a  "corner"  for  which  the  angle 
measured  in  the  conductor  is  greater  than  180  degrees.  It  is  also  evident  that  the 
current  density  is  relatively  small  along  the  left-hand  edge  near  the  centerline 
( y  =  0 )  and  near  the  top  {y  =  h).  The  small  current  density  in  these  regions 
results  from  conduction  past  corners  in  which  the  angle  is  smaller  than  180 
degrees.  We  have  discussed  how  conduction  occurs  at  these  large-  and  small- 
angle  comers  in  much  of  our  previous  work  and  we  discuss  it  in  some  greater 
detail  when  we  present  results  of  the  other  calculations. 


Isotherms  at  the  end  of  the  calculation  (i.e.,  at  2.6  ms)  are  shown  in  Figure  5.  The 
isotherms  demonstrate  the  effect  of  the  high  and  low  current  density  discussed  in 
relation  to  Figure  4.  In  particular,  there  is  evidence  of  substantial  heating  in  the 
vicinity  of  the  root  radius.  A  small  area  around  that  location  has  reached  the 
incipient  melting  temperature  of 750  K.  Similarly,  the  corners  at  the  top  and 
along  the  centerline  are  relatively  cool. 
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Figure  5.  Isotherms  in  Armature  at  2.6  ms  for  Case  1. 


Results  for  the  remaining  six  calculations  are  qualitatively  similar  and  are  not 
discussed  in  detail.  Consequently,  we  plot  for  each  of  the  calculations  only  the 
temperature  along  the  left-hand  edge  of  the  armature,  Tj,  as  a  function  of  the 
coordinate  y,  and  the  temperature  as  a  function  of  the  horizontal  distance  Ax 
from  the  edge  of  the  armature  into  the  armature  along  the  line  y  =  y^.  The  plots 
are  presented  in  groups  containing  three  cases:  Cases  1, 2,  and  3  (rg  varies);  Cases 
1, 4,  and  5  {a  varies);  and  Cases  1,  6,  and  7  (/3  varies). 

In  Figures  6  and  7  are  presented  the  plots  just  described  for  Cases  1,  2,  and  3.  As 
expected,  the  heating  in  the  vicinity  of  the  root  radius  increases  with  decreasing 
values  of  rg.  This  behavior  arises  because  as  current  is  conducted  around  a 
convex  sxurface  such  as  the  root  radius,  diffusion  arising  at  one  localized  point  on 
the  surface  does  not  overlap  with  that  arising  at  another.  Consequently,  diffusion 
away  from  the  surface  is  difficult,  the  streamlines  remain  closely  confined,  and 
the  current  density  is  high.  This  situation  may  be  compared  with  one  in  which 
current  is  conducted  along  a  concave  surface.  Then  the  skin  effects  from  localized 
points  do  overlap,  and  the  fields  and  current  diffuse  easily.  This  description  is 
analogous  to  conduction  aroimd  corners  in  which  the  angle  in  the  conductor  is 
greater  than  180  degrees  (convex  surface)  or  less  than  180  degrees  (concave 
surface),  as  explained  in  our  earlier  work.  It  is  noteworthy  from  Figure  7  that  the 
temperature  drops  to  near  ambient  at  a  distance  of  about  15  mm  from  the  edge  of 
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the  armature  in  each  of  the  three  cases.  In  addition,  the  temperature  is  nearly  the 
same  for  each  case  at  distances  of  about  5  mm  or  greater.  Clearly,  at  these 
distances,  the  details  of  the  root  radius  cannot  be  "seen";  rather,  the  current 
appears  to  have  come  from  two  lines  Lj  and  that  intersect  at  an  angle  given  by 
a+p.  The  significant  variations  in  slope  of  the  curve  for  Case  2  near  the  surface 
(Ax  =  0)  appear  to  result  from  absorption  of  the  heat  of  fusion  that  occurs 
predominantly  in  the  temperature  range  between  750  K  and  908  K,  followed  by 
the  more  rapid  heating  that  occurs  after  908  K.  The  more  rapid  heating  occurs 
because  the  resistivity  of  the  liquid  is  higher  than  that  of  the  solid.  In  Cases  1  and 
3,  the  temperature  does  not  exceed  908  K. 


Figure  6.  Final  Temperature  Along  the  Left-Hand  Surface  of  Armature  Versus  y 
for  Cases  1, 2,  and  3. 


In  Figures  8  and  9  are  shown  results  for  Cases  4  and  5,  which  correspond  to  those 
just  discussed.  In  Cases  4  and  5,  the  angle  a  was  varied  from  its  value  in  Case  1, 
while  and  /3  were  held  constant.  Clearly,  there  is  virtually  no  dependence  on  a 
for  these  cases.  Presumably,  the  radius  is  sufficiently  large  and  the  points 
and  (x^y^)  (where  the  straight  line  segments  begin)  are  sufficiently  far 
removed  from  the  line  at  y  =  Jg  that  the  changes  in  the  surface  structure  at  those 
points  do  not  affect  the  results.  Such  a  conclusion  might  not  hold  for  values  of 
sufficiently  small  that  the  origins  of  the  straight  line  segments  were  not  so  far 
removed.  Identical  conclusions  can  be  drawn  for  Cases  6  and  7  in  which  the 
angle  /3  was  varied.  Those  results  are  evident  in  Figures  10  and  11. 
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T(K) 


Figure  7.  Final  Temperature  Versus  Distance  From  Armature  Surface  Along  the 
Florizontal  Line  y  =  yo  for  Cases  1, 7.,  and  3. 


Figure  8.  Final  Temperature  Along  the  Left-Hand  Surface  of  Armature  Versus  y 
for  Cases  1, 4,  and  5. 
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Figtire  9.  Final  Temperature  Versus  Distance  From  Armature  Surface  Along  the 
Horizontal  Line  y  =  i/o  for  Cases  1, 4,  and  5. 
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Figure  10.  Final  Temperature  Along  the  Left-Hand  Surface  of  Armature  Versus  y 
for  Cases  1, 6,  and  7. 


900 


Case1,p  =  l2° 


Figure  11.  Final  Temperature  Versus  Distance  From  Armature  Surface  Along  the 
Horizontal  Line  y  =  ]/„  for  Cases  1, 6,  and  7. 


4.  Summaiy  and  Conclusions 


A  simplified  version  of  a  time-  and  position-dependent  solution  for  current  and 
heat  transport  in  a  solid  armature  has  been  solved.  Uivlike  in  previous  work,  the 
motion  of  the  armature  and  its  resultant  influence  on  current  flow  at  the  rail- 
armature  interface  were  neglected.  This  simplification  dramatically  reduces  the 
solution  time  of  the  problem  and  is  not  expected  to  affect  the  results  in  the  region 
of  interest,  namely,  at  points  removed  from  the  rail-armature  interface. 

Solutions  predicting  the  temperature  in  the  armature  at  the  time  of  exit  indicate 
that  the  angle  of  the  rear  surfaces  (within  the  rather  large  allowable  design  space 
of  8  to  24  degrees)  does  not  appreciably  affect  the  heating  at  the  root  radius.  That 
such  a  dependence  might  exist  was  a  previously  identified  source  of  concern. 
However,  as  expected,  the  generation  of  the  thermal  load  is  rather  sensitive  to 
the  radius  of  the  root  region.  The  load  increases  dramatically  for  radii  of  less  than 
about  2  mm. 

Future  work  wiU  employ  the  full  2-D  solution  for  a  moving  armature  and  rail 
conductor  system.  In  addition,  a  current  pulse  that  is  representative  of  a  more 
realistic  pulsed  power  supply  will  be  used. 
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